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Accurate analytic approximations are developed for the band gap boundaries and surface waves of 
a ID photonic crystal, making use of the semiclassical theory recently developed by the authors. 
These analytic results provide useful insight on systematics of surface states. 



42.70.Qs, 78.67.-n, 03.65.Sq, 
I. INTRODUCTION 

Aside from their intrinsic interest, surface electromag- 
netic waves (SEW) have recently been proposed [1-4] as 
a way to efficiently inject light into a photonic crystal 
waveguide, or to extract a focussed beam from a chan- 
nel. Absent some such mechanism, the insertion and ex- 
traction is problematic. In the cited works, the photonic 
crystal (PC) was a two dimensional array of rods, of infi- 
nite length normal to the plane. In this paper we consider 
SEW on a one-dimensional (ID) PC, for which we re- 
cently developed accurate semi-classical approximations. 

Surface electromagnetic waves on a ID PC were ob- 
served almost 30 years ago [5] [6]. The basic theory 
was developed at that time, [7], and can be found in 
the monograph of Yariv and Yeh [8] . More recently the 
effect of varying the thickness of the termination layer 
has been measured experimentally [9,10] and a sensor 
based on the properties of SEW's has been proposed and 
demonstrated [11]. In parallel, numerical calculations for 
SEW's in the bandgaps have been performed, and further 
aspects of the theory have also been developed [12-15]. 

Here we study the properties of semi-infinite ID pho- 
tonic crystals with termination layers of various thick- 
nesses. The dispersion equation for the SEW's is well 
known [8]. However, exact results can be obtained only 
numerically, and to our knowledge, no simple analytic 
approximations for them have been developed. Accurate 
approximations not only assist in finding exact solutions, 
but also clarify the role of the various parameters defining 
the crystal. 

Starting from the dispersion equation, we first derive 
exact expressions for the critical thicknesses at which 
SEW solutions appear and disappear. We then intro- 
duce approximate analytic forms for the dispersion rela- 
tion which are valid in specified cases. We also apply the 
semiclassical method introduced by us in [16] and [17] 
to SEW's. These semiclassical approximations provide 
simple and accurate expressions for the bandgap edges. 
They also suggest a useful empirical parametrization that 
underlies our analytical approximate forms. 

For brevity we will discuss only the case of TE waves. 



Because our methods are formally quite different from 
those recently presented in [15], in Section II we provide 
a short summary of the transfer matrix approach, in the 
notation of our previous work [16]. In Section III we red- 
erive the exact equation for SEW and obtain from it var- 
ious analytic approximations for a semi-infinite crystal. 
With these in hand, we discuss systematics of SEW's. In 
Section IV we apply the semiclassical approximations of 
[16] to surface waves, and show that the second approx- 
imation is very accurate both for the dispersion relation 
and the bandgap boundaries. This lends support to the 
parametrization introduced in Section III. In Appendix A 
we derive some closed analytic expressions for quantities 
introduced in [16] as infinite series. These allow a better 
understanding of the role of the various parameters of 
the PC. 



II. TRANSFER MATRIX APPROACH FOR A 
PERIODIC CRYSTAL 

We study surface waves of the Tamm type, that form 
at the interface between a uniform medium of low refrac- 
tive index, no, and a semi-infinite 1-D photonic crystal, 
capped by a termination layer of the same material but 
selected width. To clearly separate the periodic bulk from 
the remaining surface layer, wc split the termination layer 
of index of refraction n\ and width d c into two pieces, of 
lengths d s + d t — d c . Then the periodic array that makes 
the bulk of the crystal consists of "cells" each made of 
three uniform layers of widths dt, g?2 and d\ — d t whose 
respective refraction indices are n\ , and n\. The 
cells are reflection symmetric when 2dt = d\ . A sketch is 
shown in Fig. 1. The initial cell, extending from z = 
to z = d = d± + c?2 will be given the index 0, the second 
index 1, and so on, so that the p-th cell extends from pd 
to (p + l)d and has n{z) = m when pd < z < pd + d t 
or pd + d t + g?2 < z < (p + l)d and n(z) = n-i when 
pd + d t < z < pd + d t + g?2- As is customary, we will 
suppose ?ii > ri2 > rig. The rest of the cap layer extends 
from z = —d s to 0, and the uniform medium is located 
to the left of z = —d s . 
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We choose a coordinate system in which the layers 
have normal vector along OZ. An obliquely incident plane 
wave defines the OX axis. For monochromatic TE waves 
the electric field is parallel to the OY axis. As in [16], we 
write 

V = E{z)e y e l ( k l 3x - ut ) 

H=[H x (z)e x + H z (z)e z ]e^ x -^) , (1) 

where lu is the angular frequency, k = ui/c is the vacuum 
wavenumber and (3k is the (constant) x-component of 
the wavevector of modulus k{z) = um(z)/c. For a TE 
wave entering the 1-D photonic crystal from a uniform 
medium, one has 



(3 = tiq sin 6*o , 



(2) 



where 8 is the angle of incidence measured from the 
normal. Maxwell's equations require that in a periodic 
medium the E(z) introduced in cq. 1 satisfies 



d 2 E(z) 
dz 2 



+ k 2 (n 2 (z) - f3 2 )E(z) = 



(3) 



The solutions are well known. In the p-th cell, the electric 
field corresponding to TE waves can be written as 



E{z) = a p e lkliz - pd) +b p e- lkl 



(z-pd) 



(4) 



when pd < z < pd + d t , and fc x = k^n\ — (3 2 = knip. 
Similarly 



E(z) = c p e lk2< - z - pd ^ + d p e- lk2i - z - 



P d) 



(5) 



when pd + d t < z < pd + d t + d 2 , and k 2 = kyjn 2 , — (3 2 
kn 2 p. Also, 



E(z) = e p e 



— P Jki(z-pd) 



+ f p e 



-ik\ (z—pd) 



(6) 



when pd+d t +d 2 < z < (p+l)d. Matching these solutions 
and their derivatives at the two interfaces, one finds the 
transfer matrix, M. [beware: this matrix is called M.~ x 
by some authors]. 



& p+i 



M 1 bp 



with 



A = e lkldl (cosk 2 d 



A B 

B* A* 



^K£ + l )si,)/, -" /2 



P ) , (7) 



B = e^- 2 ^U^~^) S mk 2 d 2 . (8) 
2 \ki k 2 J 

The periodic (Bloch) waves of the infinite crystal are 
eigcnsolutions of the transfer matrix 



( a P+i \ = p ±i<t> ( a v\ = 
\b p +i 



(9) 



and therefore the secular equation is 

(A - e ±i<t '){A* - e ±i< ^) - \B\ 2 = 
with eigenvalues 



= Re(A) ± ^Rc(A) 2 - 1 , 



(10) 



(11) 



where we have made use of det M. = \A\ 2 — \B\ 2 = 1. 
The corresponding eigenvectors are, up to a normaliza- 
tion factor 



B 

,±i<t> _ A 



(12) 



In allowed bands <f> is real and the bandgap boundaries 
are defined by the condition RcA = ±1. In bandgaps 
the eigenvalues X± = e ±l ^ are real since \RcA\ > 1, and 
therefore 



A± = Rc(A) ± v/Rc(A) 2 - 1 



(13) 



with A_A + = 1. We choose the solution that gives a 
damped wave when z — > +oo. It is A_ ( A + ) when 
Re(A) > 1 ( < — 1.) We write it as simply A and bear in 
mind that |A| < 1. 



III. SURFACE WAVES 

These are states which decay in both directions, as one 
moves away from the surface of the photonic crystal. To 
the left of z = we have the variable portion of the cap 
layer, and the uniform external medium. The electric 
field is written as 



E(z) = a s e lklZ + b s e- tklZ 
when — d s < z < and 

E{z) = b v e q ° z , 



(14) 



(15) 



with g = k\J fi 2 — ng, when z < —d s . Matching at the 
interfaces, and choosing the damped wave solution for 
z > 0, one finds 



ki 



X-A-B 



X-A + B 

where for simplicity we absorb a phase into 

B = e - 2lklds B . 



(16) 



(17) 



Eq. 16 determines the dispersion relation k = k{(3) for 
the surface waves. It has to be solved numerically, and 
we will refer to the solutions thereby obtained as "exact" . 

We begin by examining the structure of eq. 16. From 
the definitions of go and k\ given above, we see that the 
left hand side depends on (3, but not on k. Eq. 10 shows 
that |A — A\ = \B\ = \B\, so that by writing 
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@\—A = arg(A - A) , 9^ = arg(S) , (18) 
eq. 16 becomes: 

o i0\-A _ gi#B 



9o = _f 



tan 



'\—A 



fci _ ' e^A-^i + e i8s y 2 

Next we look at arg(B). From eq. 8 we note that 

_ ki _ n2f3_ _ nip_ < 
h k 2 nip n 2 p 

according to the choice n\ > n 2 made earlier. Hence, 



(19) 



(20) 



7T 



dg = — + kn 1(i (d Y - 2d c ) + ■ 



(21) 



with 4> s chosen to be or —ir depending on the sign of 
sin k 2 d 2 . As shown in [15], this sign is characteristic of 
each bandgap, unless the latter shrinks to zero width in 
what is called an optical hole [15]. Defining 



9(0) = tan" 



we can rewrite eq. 19 as 




(22) 



0\- A (k) = 29(0) + knx^dx - 2d c ) + (j) s + - + 2iw 

(23) 

with v = 0, ±1,±2, The l.h.s. and the second term 

on the r.h.s. of this equation depend on k, while the 
others do not. 

The term 2nv arises from the inverse of the tangent 
appearing in eq. 19. When for a given and d c = d c , 
the v = solution is k = ko, one sees easily that k = 
ko is also a solution corresponding to the same and 
d c = d c fi + vif / (n-\_f) ko) i with v = ±1,±2, • • •. This is 
analogous to the well known property of solutions of the 
Schrodinger equation for a finite square well potential: 
increasing the width of the well by a half-wavelength, 
produces a state of the same energy but one additional 
node. For simplicity we will focus, from here on, on the 
case v = 0. 

If 6\—A was linear or quadratic in k for fixed 0, we 
could easily solve eq. 23 for k(0, • • •) in terms of the other 
parameters, and so identify the SEW's. Since that is not 
the case, one is reduced to numerical or graphical meth- 
ods of solution. To see how this works, we consider an 
example in the first bandgap, taken from [16]: a PC with 
parameters n\ — 2, d\ = 100 nm, n 2 = 1.5 and d 2 = 250 
nm. In Fig. 2 we plot separately the left hand side (con- 
tinuous line), and the right hand side for several values 
of d c (dashed lines). In the first bandgap, sinfc2C?2 > 
and therefore we have set <p s — — n. As k varies over the 
bandgap from kh to ka, the argument 6\-a increases 
from — 7r/2to7r/2. This is a generic feature, as discussed 
in [18]. The intersection of any dashed line with the con- 
tinuous line defines a corresponding solution for k. One 



sees that as d c decreases the corresponding k increases, 
as expected (think of the analogy with the solutions of 
the Schrodinger equation). The extreme values of d c for 
which a solution exists will therefore be those for which 
the r.h.s. of eq. 23 takes the values ±7r/2. This gives the 
values 



dr. 



di 



kRn h p 

1 



2 k L ni,p 



e(/3). 



(24) 



In the example shown in the figure, corresponding to /3 = 
1.2, one finds d Cym i n = —27. nm and d c . max = 86.5 nm. 
The negative sign merely indicates that there is a surface 
wave solution for d c ranging from to d c , max . 

To proceed further in the analysis of the SEW solutions 
requires values for the band edges fci, fcjj. These will 
be obtained from our semiclassical approximation [16] in 
Section IV. Before we delve into that, we first introduce 
an empirical approximation to 9\-a which will be justi- 
fied by the semiclassical theory. For the first bandgap, 
we write 



(9 



(e) 
X—A 



sin 



w/2 



(25) 



with k m = (kfj + fcf,)/2 and w — k^ — k^. In Fig. 3 we 
compare the exact and the empirical forms of 9\-a for 
four values of j3 ranging from no to n 2 - It can be seen 
that the approximation is quite satisfactory. 

Based on the above development, we now derive two 
analytical approximations for the dispersion relation. 
The first is valid when the crossing point in Fig. 2 lies in 
the linear portion of 0\-a, while the second is valid near 
the upper and lower ends of the curve. 



A. Solutions in the middle of the bandgap. 

These are of particular interest because the damping is 
stronger, strongly confining the wave to the surface [9]. 

Assuming that k — k m « w we can approximate 
sin _1 (2(fc — k m )/w) ~ 2(k — k m )/w and eq. 23 then 
has a solution 

k „ &{f3) + k m /w + 4> s /2 + K{v+l/A) 
l/w + n w {d c -d 1 /2) 

In this way we can easily construct k(/3) for a fixed value 
of d c , or conversely, study k — k(d c ) for fixed (3. The role 
of the bandgap parameters k m and w is also quite easy to 
analyze. Fig. 4 shows the accuracy of this approximation 
when d c = 25 nm. For this example, when /3 > 1.4 the 
approximation ceases to be valid and one has to resort to 
other approximations described in the next subsection. 
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B. Solutions near the bandgap boundaries. 

These approximations will be useful in analyzing the 
results of experiments like those of Robertson and May 
[10], where the SEW's appear very close to the bound- 
aries. As seen in Fig. 3, the linear approximation to 
the arcsine fails near the band boundaries. We discuss 
solutions near the lower boundary, but similar approxi- 
mations can be developed for the upper boundary. 

When k is slightly above kh — k m — w/2, it is conve- 
nient to introduce £ > via 



Then 



sin 



k kf: 



\ ( k kf t 



-d-c). 



(27) 



w/2 



7T 
2 



2C 



7T 

-2 +2 V „■ 



k-kr 



Inserting this into eq. 23 gives 



k — k L 



+ p(k-k L ) = A 



(28) 



(29) 



with 



p = n ll3 (2d c - di) 

A = 2Q + <j) s + it{2v + 1) - k LP 



and solving for k — k^ 



k = k L 



p< 




(30) 



(31) 



which is the desired solution, k = k{(3), near the lower 
bandgap boundary. Fig. 5 shows an example of the ac- 
curacy of this expression. Furthermore, when pA is small 
compared to 1 /w one can expand and find 



k = k L + -wA 2 
4 



(32) 



have an optical hole in the second bandgap, so we exam- 
ine this situation to clarify what that entails. 

In eq. 23, the angle 6\-a ranges from — tt/2 + n(q — 
1) to 7r/2 + ir(q — 1) in the q-th bandgap. Following 
the previous argument we find the critical thicknesses for 
appearance/disappearance of SEW's to be 



di 
2 



,:,„,„ - — v -j—^- — ( 6(/3) + ^ - ^(q - 1) + 7H/ 
2 KRn\,p V 2 2 



1 



,,,„,, - - - , 0(/?) + ^- ZL(g_2)+7H/ 

2 k]_,ri\_jj \ 2 2 



(33) 



In the second bandgap, Fig. 6 shows that an optical hole 
occurs, where the width of the gap shrinks to zero. (Sec 
Section IV and [15] for the exact location.) At the optical 
hole, fc 2 c?2 = gi vm g a change of sign of sin/c 2 ^2 at this 
critical [3 — (3 h- Correspondingly, one finds that 4> s = 
for (3 < (i h and <j> a = —it for (3 > [3 Q h- The continuous 
lines in Fig. 7 show d c ^ max and d c . m i n , computed from eq. 
33 with v = 0. The horizontal dashed lines correspond to 
various values of d c . The plot shows that for d c = 10 nm, 
there is a SEW for (3 less than approximately 1.25, and 
that the SEW reappears for (3 > (3 h- Similarly, for d c = 
15 nm there is a SEW when (3 is less than approximately 
1.325, reappearing for (3 > (3 oh . For d c = 20 nm SEW 
appear for (3 < (3 h , and for (3 greater than approximately 
1.4. Finally, for d c = 40 nm there is a SEW only when 
13 < [3 oh . 

Fig. 7 shows surface wave solutions when d c = 0.15<ii 
and d c = 0.40 d\, and confirms the above discussion. 
Note that numerically it is rather difficult to locate the 
end point of the surface wave solution for d c — 0.15 d\, 
because the solutions run very close to the band bound- 
ary. In contrast, eq. 33 and the graph Fig. 7 locate that 
end point very easily. 



D. Solutions when (3 > ni 



In this regime, ki — k^Jrii, — (3 2 becomes imaginary 
and we write k-2, = iq2- The expressions for A and B 
become 



which again manifests the dependence of k— fc^ on d c and 
w, and allows one to construct k = k{(3) very easily. Fig. 
5 shows again the validity of this approximation. Note 
also that the condition A = determines the location of 
the zone boundary. Writing it out, one recovers eq. 24, 
so that eqs. 31 and 32 do not violate this exact relation. 



C. Surface states in the second bandgap 

In the examples discussed above we focussed on disper- 
sion relations for SEW in the first bandgap. The specific 
system considered has no optical holes. However, it does 



B 



* e i fcl (*-2*o (h + «i) sinhQ2d2 . 



2 V 92 ki 



(34) 



As (3 increases, the hyperbolic functions soon become 
very large, giving 



A ~ 1 f ,ikid 1 +q 2 d 2 (l, l ( kl 92 

A -2 G l 1+ 2U k[ 

i e ik 1 (d 1 -2d c )+q 2 d 2 + 92 

4 V 92 h 



B 



(35) 
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and therefore: 



with 



Re(A) ~ -e« 2d2 (cosfcidi - Tsinftidi) 



1 Ai <? 2 \ _ 1 + - 2/3 2 



2 V?2 fci/ 2 ^/(nf-/? 2 )^ 2 -^) 



(36) 



(37) 



which is independent of fc. The stationary points of the 
quantity in parentheses on the right of eq. 36, are at 
tanfcidi = — r. At these points: 



Re (A) = (-)«ie 92 *Vl + r 2 



(38) 



Therefore Re(A) alternates in sign from one bandgap to 
the next, and the amplitude of the oscillations is large in 
most of the range of values of q 2 due to the exponential 
factor. Writing Rc(A) = (-l) q \Re(A)\, and using cq. 13, 
we obtain 



A = (-)« (\MA)\-V(MA)) 2 -l) 



(-)« 



2\Re(A)\ 



, (39) 



which is much smaller than \A\ over most of the bandgap. 
Therefore 



y\-A 



~ arg(— A) = — 7r + k\di + tan 1 T , 



(40) 



where the last term is independent of k and therefore 
9\-a becomes linear in k. Figure 9 shows the accuracy 
of eq. 40: beyond (3 — 1.6, the exact (continuous) and 
the approximate (dashed) dispersion relations practically 
coincide. Using this approximation and equating 6\-a 
to (q — l)n ± 7r/2 we get an expression for the bandgap 
boundaries 



nipdi 



an ± tan 1 T 

2 



(41) 



Figure 10 shows that the band boundaries are quite well 
reproduced by this approximation. The exception is the 
lowest boundary when (3 < 1.6. We can also use the ap- 
proximation for Q\—a to predict the dispersion relation: 
inserting 40 into eq. 23 and solving for k, one finds 



k 



1 



ll/3«c 



e 



7T (t>e 1 1 „ 

3 4 + y-2 tan r 



(42) 



As figure 10 shows, the exact and the approximate curves 
for k = k(fi) are very close (indistinguishable when (3 > 
1.6). Again, eq. 42 shows very explicitly the role of d c 
and of the indices of refraction in determining k. 



IV. SEMICLASSICAL APPROXIMATIONS 



A B 

B* A* 



1 

tb 



(43) 



that is, A = l/t* B and B = —Ar* B . The semiclassical 
expressions for r B and t B are found in eqs. 35 and 36 of 
[16] for the first and second approximations discussed in 
that reference. 



A. First approximation 

In this approximation, and for a single cell, eq. 35 of 
[16] gives 



-s* sinh(7i cf) 



71 cosh(7i<i) — iS q sinh(7id) 



t (1) - 



R 9 7i 



71 cosh(7irf) — iS q sinh(7irf) 



(44) 



Therefore 



A^ = (-) g ^cosh7ic? + z— sinh7irf 

B W = (_) 9 fi sinh7ld ) 
7i 



(45) 



where the subscript s stands for for semiclassical, and 
the superscript (1) denotes the first approximation in- 
troduced in [16]. The explicit expression for s q is given 



below, eq. 46, 71 = yj \s q \ 2 — <5 2 and the detuning from 

the q-th Bragg resonance is S q = n av j](k — k B ), with 
k B = qn / (n av jjd) . In addition, 71c? is the exponent of 
the damped wave constant, A, as we can easily confirm: 
inserting RcA^ into cq. 13 we find = (-)^e-^ d . 

Using expressions from [16], for an asymmetric cell we 
can write s q in compact form 

S q = - l - In 1 M e ii4>i+4>2)/2 sin((( ^ _ fayy (4g) 

a n 2 /s 



with (f>i = 2(— k\ + S q )d t 

(f> 2 = 2(-k 1 d t - k 2 d 2 + S q (d t + d 2 )) ■ 

From eq. 45 we find 

- - (-)(«- 1 ) sinh( 7l d) ^1 + ■ 



7i 



(47) 



(48) 



and therefore # A (i)_ A (i) 



= arg(AW - A&) = (q - 1)tt + tan" 1 ^ 

7i 

= {q- 1)tt + sin" 1 A ~ ( q - 1)^ + sin" 1 A- , (49) 



The elements of the transfer matrix can be related to 
the amplitudes for transmission ts and reflection r B , for 
a single cell: 



where in the last step we have replaced the slowly varying 
|s g (fc)| by Sb, its value at the Bragg resonance k = k B . 
Since by definition, S q = (k — fcs)n a „ ;/ 9, we arrive at 
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id)- 



(i) ~ (q - l)w + sin 1 



k — k B 
w/2 



(50) 



with w = 2S B /n av ^. This has the same form as the 
empirical parametrization used above, but provides ex- 
plicit estimates for the width and position k m — k B of 
the bandgap. 



B. Second approximation 

The second, and more accurate, approximation intro- 
duced in [16] leads to similar but more involved expres- 
sions for and B^ . From eq. 36 of that reference, 
wc find 



A?) = {-) 



cosh 72 d 



[(i + H')»- W)] Slnh72 , 

(1 - M 2 ) 7 2 



(51) 



where 72 = \J\s q \ 2 — J] 2 and r\ q = S q — ic 2 . In the ap- 
pendix we give new analytic expressions for u = V\ (z = 0) 
and for C2, that sum the series written in [16], eqs. 18 
and 19. As in the previous case, one easily finds that 



A 



(2) 



(-)'e" 



-72 d 



(2) (2) 

Inserting A s and Bs into cq. 16, and solving, one 
finds the corresponding predictions for k = k{[3). Before 
presenting these results, we will show the usefulness of 
this second approximation in giving accurate values for 
the bandgap boundaries. 

In eqs. (31) and (32) of ref. [16] we showed that the first 
approximation provides simple estimates for the bandgap 
boundaries. In the second approximation, the condition 
72 = defines the band boundaries, because it corre- 
sponds to infinite decay length of the surface state. Us- 
ing the explicit form of 72 from eq. (25) of [16], we can 
write this condition, for the g-th bandgap boundaries, as 



~\s q (k L )\ = (k L - k B )n av ^ + r 2 
\s q {k R )\ = {k R - k B )n av ,p + r 2 , 



(52) 



where T2 = — ic 2 . The dependence of |s g | on k is fairly 
smooth as long as k remains within a bandgap. To a 
good approximation one can expand around k — k B and 
write 

\s q (k)\ =S B + £(fc - k B ) - ^ V (k - k B f , (53) 



with 



S B = \s q (k B )\ = ~An ( — ] SllldiJ . 

a \^2, 



a B = irq 



P 



n2,fjd 2 
n av ,f3d 
did 2 



(n 2 ,p - ni i|S ) 

1, ( n hp\ 
: In I I p cos as 



d \n 2 



1 , / ni,0 \ 
11= — In p sinas 

d V n 2,/3 



(54) 



Inserting this expansion into each line of 52 leads to the 
desired analytic expressions 



k L = k B + - 
V 



k R = k B + 



- yj {n a v,(i + 2 + MSb + r 2 ) 

- n av ,p + £, 



with r 2 evaluated at k = k B : 

di - d 2 ( . ( nx.fi 
r 2 = — r^— I In ' 



Ad 2 



112,0 



sin(2a_e) 



(55) 



(56) 



In Fig. 6 we compare the exact bandgap boundaries to 
those of eq. 55 and, as can be seen, the latter works very 
well except very close to [3 — n 2 , where the semiclassical 
approximation is expected to fail. For (3 close to n the 
values of 77 are small and one can neglect them in cq. 53. 
This leads to simpler forms: 



k T , = kt 



kn = k B + 



S B + r 2 

n a v,p + £ 
S B - r 2 
n av ,p - £, 



(57) 



Neglecting the contribution from r2 we recover the ex- 
pressions of the first approximation (71 = 0) already dis- 
cussed in [16]. 

Direct inspection of the expression for Re(A) derived 
from eq. 8 shows that the points where the bandgaps 
shrink to zero width correspond to special values of the 
optical depths of the layers: k\d\ — miir, k 2 d 2 — m 2 ir, 
with mi, m2 integers. These give 



n lt/ 3d 1 m 2 = n 2 ^d 2 mi , 



so that: 



n\d\m2 — n\d\m\ 



d\m\ 



(58) 



(59) 



and k a h — miTr/(dini.fj oh ) 1 in agreement with [15]. The 
point where the second bandgap shrinks to zero width, in 
Fig. 6, corresponds to mi = m 2 = 1, q = 2. It is easy to 
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check that with the above semiclassical approximations 
one finds exactly the same optical hole. Note that at this 
critical value 



a B = nq- 



k 2 d 2 



irq- 



m 2 



= 7T 



(60) 



(61) 



kidi + k 2 d 2 'mi+m 2 
so that r 2 = 0, Sb = V = 0, and therefore 

k L = k R = k B ■ 
C. Results for SEW's 



In Fig. 11 we compare the exact and semiclassical re- 
sults for k = k(/3), choosing three thicknesses for the cap 
layer: d c — dt = 25, 50 and 75 nm. The first approxima- 
tion becomes inaccurate when (3 exceeds 1.4, but gives 
accurate results up to that value. The second approxi- 
mation is so close to the exact values that one can see the 
difference only for values of (3 very close to the critical 
value (3 = n 2 = 1.5. Beyond that, our semiclassical ap- 
proximations cannot be applied, since k 2 becomes imag- 
inary. In the second bandgap, the accuracy of the first 
approximation is significantly worse, whereas the second 
approximation is as good as for the first bandgap. For 
brevity, we do not show any figures for A and B as func- 
tions of k and (3. For most of the values in the first 
bandgap the agreement is similar to that seen in fig. 11 
for the dispersion relations. 



V. SUMMARY AND CONCLUSIONS 

By considering the semi-infinite limit of a ID PC 
we have derived a dispersion equation for SEW's, valid 
for termination layers of selected width. Our goal has 
been to clarify the systematics of the solutions of this 
equation. To do so we first discussed a graphical so- 
lution. This allowed us to derive analytic expressions 
for the critical thicknesses at which solutions appear at 
the bandgap boundaries. Further, by introducing a suit- 
able parametrization of 0\-a, 25, we have derived simple 
approximations for solutions cither in the middle of the 
bandgap or near the edges: eqs. 26, 31, and 32. We 
tested them by an example whose first bandgap has no 
optical holes. We then extended the analysis of criti- 
cal thicknesses to the second bandgap, where there is an 
optical hole. The appearance and disappearance of solu- 
tions as a function of cap layer thickness is again easily 
predicted: eq. 33. For completeness we have also ex- 
amined solutions when (3 > n 2 , and again found very 
simple approximations eqs. 41 and 42, valid over most 
of the range n 2 < (3 < m, that allow us to study the 
systematics in a transparent way. 

Finally, in section IV, we applied the semiclassical ap- 
proximations derived in [16] to SEW's. The second ap- 
proximation is very accurate in predicting the bandgap 



boundaries and the dispersion relation. We took advan- 
tage of this to derive accurate and simple approximations 
for the boundaries: eqs. 56 and 57. The first approx- 
imation already supports the validity of the empirical 
parametrization of 9\^a- cf. eqs. 25 and 50. 

In conclusion: we have presented a set of analytic re- 
sults, exact and approximate, that clarify the systemat- 
ics of solutions for surface EM waves in semi-infinite ID 
photonic crystals. In addition, we have found simple ana- 
lytic expressions for the bandgap boundaries that should 
be useful for the design of PC configurations. We plan 
next to extend our results to layered configurations with 
cylindrical symmetry [20,21]. It would also be interest- 
ing to see whether the systematics found here apply to 
surface states in 2D and/or 3D photonic crystals. 
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APPENDIX A: CLOSED FORMS FOR SERIES 
USED IN THE SEMICLASSICAL 
APPROXIMATIONS 

1. Analytic expression for u = vi(0) 

The function v\(z) is defined in cq. 19 of [16] as 

i2ir(m—q)z/d 



id x 
9^ 2-^, 



27r ^— ' m—q — 5 a d/ / n 



(Al) 



with s m given in eq. 12 of [16]. In the case of interest 
here, the refractive index is piece-wise constant, so that 

_ J_ \p i n ei z 3 + °) 2i{-i 1 {z j )+kn av ^z j --Kmz j /d) 

Sm ~ 2d ^ n e ( Zj - 0) 



(A2) 



and n e (z) = \J n 2 (z) — f3 2 . Therefore 
Vl (z) = - V e 2i ( - ^') + kn av,0Zj - irqz 3 )/d 

47T ^ 

i 

n e {z d + 0) 



In 



e i2TT(m-q)(z-Zj)/d 
E ■ (A3) 



n e (zj — 0) *—f m — q — S a d/7t 



Using eqs. 1.445.7 and 8 of [19], we find 

Vl {z) 



4^~n e (^-0) 



S q d 



e i[(2(z-Zj)/d-l)S q d] 

sin S q d 



(A4) 
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with z = z when < z — Zj < d and z = z + d when 
z — Zj < 0. 

For the asymmetric cell of the photonic crystal of interest 
and 2 = 0+, this gives: 



u = ui(0+) 



— e 



4 n lf} 
1 

Sqd 



1 



,i(l-2d t /d)S q d 



S q d sm(S q d) 



J(l-2(d 2 +d t )/d)S q d 



sm(S q d) 



with </>i and </>2 given in eq. 47. 

2. Analytic expression for C2 

According to eq. 18 of [16]: 
id 



9.tt ^— ' m. — a — 



2ir ^— ' m — o — 5 a d/ir 

m=tq H 

and using eqs. 46 and 47, we have 

2 



c 2 = Afln'^ 



sin 2 (P™ + T) 



So with 



m^q 



m — q — e 



where 



(A5) 



(A6) 



(A7) 



P = ^-^ . r=— i-i(n2,/j-ni,^) (A8) 

and e = Sqd/ir. Using eqs. 1.445.6 to 8 of [19] we arrive 
at: 
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have the desired analytic expression for c 2 . 
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FIG. I. Sketch of a typical ID PC, for the case used in calculations (rn = 2.0, di = 100 nm, 712 = 1.5, d,2 = 250 nm,) A 
symmetric unit cell was chosen, and d c = 80 nm. 
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FIG. 3. Exact (continuous line) v.s. empirical approximation (dashed line) to the argument of A — A. From left to right, 
P = 1.1, 1.2, 1.3 and 1.4 
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FIG. 4. Dispersion relation k = k{0) when d c — 0.25di. Continuous line: exact solution of eq. 16. Dotted line: linear 
approximation, eq. 26, based on the empirical form of arg(A — A). The dashed lines are the boundaries of the first bandgap. 
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FIG. 5. Dispersion relation when d c — 0.75di . Exact (continuous line) v.s. approximations based on the empirical form of 
arg(A — A). Triangles: eq. 32; filled circles: eq. 31. The dashed line is the lower boundary of the first bandgap. 
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FIG. 6. First and second bandgap boundaries: Thin continuous lines: exact. Dashed lines: second semiclassical approxima- 
tion described in text. 
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FIG. 7. Second bandgap: Continuous lines d^ min and d c , max as predicted from eq. 33. Dashed horizontal lines: from bottom 
to top d c = 10, 15, 20 and 40 nm. 
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FIG. 8. Second bandgap: Continuous lines: zone boundaries. Dispersion relation for d c = 0.15di (dashed line), and 
d c = 0.40di (dotted). 



12 



0.01 0.02 0.03 0.04 0.05 0.06 

k[ nm- 1 ] 

FIG. 9. Argument of A — A in first bandgap. From left to right: /3 = 1.51, 1.6, 1.7, 1.8 and 1.9. Continuous lines: exact; 
dashed lines: linear approximation of eq. 42 
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FIG. 10. First and second bandgaps when (3 > ni. Continuous lines: exact boundaries. Dashed lines: approximation 
described in text. Surface waves: solid line with + symbols: exact; dotted lines: approximation given in the text. 
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